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We study the dynamics of states perturbatively expanded about a harmonic system of loop quan- 
tum cosmology, exhibiting a bounce. In particular, the evolution equations for the first and second 
order moments of the system are analyzed. These moments back-react on the trajectories of the 
. expectation values of the state and hence alter the energy density at the bounce. This analysis is 

' performed for isotropic loop quantum cosmology coupled to a scalar field with a small but non-zero 

constant potential, hence in a regime in which the kinetic energy of matter dominates. Analytic re- 
strictions on the existence of dynamical coherent states and the meaning of semi-classicality within 
these systems are discussed. A numerical investigation of the trajectories of states that remain 
semi-classical across the bounce demonstrates that, at least for such states, the bounce persists and 
, that its properties are similar to the standard case, in which the moments of the states are entirely 

' neglected. However the bounce density does change, implying that a quantum bounce may not be 

guaranteed to happen when the potential is no longer negligible. 
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I. INTRODUCTION 



Loop quantum cosmology [l| imports mathematical constructions and results obtained in the general setting of 
quantum geometry realized in loop quantum gravity into quantum cosmology. As a result, the quantum repre- 
^ ' sentation is unitarily inequivalent to the one used in Wheeler-DeWitt quantizations of cosmological models, a formal 
difference entailing also dynamical changes in the behavior of models for the very early universe. Loop quantum cos- 
^ mology thus allows one to explore the physical implications of quantum geometry and its unconventional space-time 
Q\ ■ structures in a tractable setting. 

CO ' One of the main changes implied by loop quantum gravity is the use of holonomies , objects of the form exp(icurv) 
. replacing the usual curvature expressions "curv" (extrinsic and intrinsic curvature, or space-time curvature compo- 
' nents) in the classical equations. The form of the curvature functional "curv" appearing in holonomies depends on the 
' specific quantization of the Hamiltonian constraint used [6^ as well as on details of the reduction to a symmetric con- 
text within quantum gravity [7]; however, within isotropic cosmological models it is a function only of the scale factor 
and its time derivative. Moreover, neither is the full theory of loop quantum gravity uniquely and unambiguously 
formulated yet, nor is the reduction from general physical states to isotropic ones fully understood at a quantita- 
tive level. Despite such ambiguities, several general conclusions can be inferred from mathematical properties of the 
curvature-holonomy replacement, most importantly related to the boundedness of holonomy expressions in contrast 
to the curvature components. For instance, the Friedmann equation relates extrinsic curvature act of spatial slices in 
a dynamical isotropic universe to the energy density of matter. If extrinsic curvature in the equation is replaced by a 
bounded function of a certain form, boundedness of the energy density automatically follows in spite of its classical 
divergence at the big bang singularity. 

Replacing curvature expressions by bounded functions is not simply done by hand but is motivated by quantum 
theory. Quantization also means that one is dealing with states rather than classical geometries subject to the 
usual evolution equations. Fundamentally, the use of holonomies primarily implies that the dynamics of quantum 
states of universe models is governed not by the differential Wheeler-DeWitt equation d, |9|, but by a difference 
equation [l^, [ll| . These discrete dynamics turn out to be non-singular [H, : wave functions extend uniquely across 
the classically singular big bang. But while the use of holonomies and the discreteness of the difference equation 
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implied by it are important for this result, detailed properties of the transition depend on various kinds of quantum 
effects that the use of the wave function introduces. In particular, in addition to quantum geometry effects such as 
discreteness, there are in general quantum back-reaction effects that characterize deviations of the quantum evolution 
of expectation values from the classical evolution caused by a changing wave function. Simple intuitive boundedness 
results on energy densities can be established only if quantum geometry effects can be shown not to be overpowered by 
potentially adverse quantum back-reaction effects. In loop quantum cosmology, this competition of different quantum 
effects has not been explored sufficiently in generic regimes. (Examples in Wheeler-DeWitt quantum cosmology have 
shown the sensitivity of some quantum effects related to singularity avoidance to state properties [14 16J.) 

Typically attention has been focused on harmonic models, in which the quantum back-reaction is entirely absent. 
This has the advantage of being analytically tractable and has been shown to generically lead to bouncing solutions. 
The robustness of these results for models that deviate from these special cases has only begun to be explored and 
in this paper we further develop a systematic perturbation theory approach to this question [TtI . [18| . Key results 
reported here include the interplay of reality conditions (physical normalization of states) and uncertainty relations, 
which in non-harmonic cases turns out to be important for consistent evolution, numerical or otherwise. By numerical 
studies, we demonstrate the sensitivity of bounce properties such as the density to quantum effects. Also the high 
sensitivity of evolution thr oug h the bounce on initial values, observed for harmonic models based on an analysis of 
dynamical coherent states (T9[[20|. is qualitatively confirmed in our more general setting. Although the leading-order 
perturbations used here do not allow quantum dynamical corrections large enough to (potentially) remove the bounce 
induced by the holonomy modification, the sensitivity we find appears sufficiently strong to warrant caution about 
extrapolations of the kinetic-dominated behavior to models with significant matter potentials. 



According to the holonomy replacement scheme, one of the basic objects in loop quantum cosmology is the holonomy 
exp(zcurv) with an expression "curv" linear in curvature components. For an isotropic cosmological model, this means 
curv = g{a)a in terms of the scale factor a, with an a-priori arbitrary function g{a). Since holonomy operators in 
loop quantum gravity are related to the discrete quantum geometry, the specific form of g{a) depends on how those 
structures are realized and how they reduce to the case of isotropic models, which is referred to as lattice refinement 
[21II22II. At present, g{a) does not appear to be fixed uniquely by reduction schemes from the full dynamics, and so we 
will parametrize it as a power law, g{a) — goa^^ with a dimensionless number x and a constant 170 whose dimensions 
are such that g{a)a in the holonomy exponential is dimensionless. The dimension of go thus depends on the value of x. 
(Several types of phenonienological analysis, mainly in the context of early-universe cosmology, have put restrictions 
on the ranges of parameters [23l - t27| ) . 




Holonomies are used in loop quantum gravity not at the level of equations of motion but at the level of constraints 
which determine the dynamics of phase-space points. In the present context (spatially fiat isotropic models), we have 
canonical variables c = 7a and its conjugate \p\ = a?: {c,p} = 87r7G/3Vo with the Barbero-Immirzi parameter 7 
(28I [29| , the gravitational constant G and the coordinate size of a finite region chosen to integrate the full symplectic 
structure specialized to homogeneous fields. (Geometrically, the momentum p actually refers to a densitized triad 
and thus can take both signs, depending on the orientation of the triad. In this article, we will assume p to be 
positive without loss of generality.) The parameter Vq is made implicit by absorbing it into the phase-space variables, 

1/3— 2/3~ 

c := Vq c and p := Vq p, but one still must ensure that final results depend only on ratios of c and p independent 
of rescaling by Vq. See Appendix iBl for an explicit discussion of this so-called fiducial volume. 

With these canonical variables, the holonomies we will use are exp(i/(p)c), where f{p) — /oP^ is obtained simply by 
rearranging the previous parameters in g{a): for fop^c = f{p)c = g{a)d = goo^^d, we identify /o = 7~^V'q ^^^^^''^^5o- 
With this Vb-dependence of /o we ensure that the product f(ji)c appearing in holonomies is independent of Vq, in 
addition to being dimensionless. It will be useful to apply a canonical transformation such that the combination f{p)c 
appears as a canonical variable, which must then be conjugate to 



II. LOOP QUANTUM COSMOLOGY 




A. Canonical variables and a harmonic system 



V := 



3p 



3„2(l-x)v2(l-x)/3 



such that {foP^c, V} = 1 . 



(1) 



8^7G(1 - x)f{p) 



87r7G(l - x)fn 



Since f{p) scales by a factor of V V cx p/ f{p) scales by a factor of Vq if Vq is changed. The volume of our fiducial 



3 



cell, vol = a'^Vo, is then obtained as: 



voi.(5!^2^ii^y)^ . (2) 

In addition to f{p)c used in holonomies, another combination independent of Vq is f{p)^/p — /op'-^^^^-'/^. The physical 
meaning of this quantity is the lattice spacing L in an underlying discrete state that gives rise to the dynamics captured 
in the isotropic model. In terms of this parameter, we can write holonomies as exjpiijL'H) with the Hubble parameter 
H. The lattice spacing can also be obtained as the ratio vo\/V = 87r7G(l — x)L/3, written as 

vol 87r7G' (1 — a;) ~ . i+2x 
^ = ^ foP ^ 

in terms of the Vb-independent /o := foVg^^^^'^^'^ . This equation shows that L is independent of the scale factor only 
for a; = — i. 

Via holonomies, the combination f{p)c appears in dynamical equations, chiefly the Hamiltonian constraint which 
follows from multiplying the Friedmann equation by . In canonical variables, this provides 

1 9 r, — r SttG „ 
- ^c^/H + — = 3 

7"^ 3 

with the matter Hamiltonian -ffmattcr typically depending on p as well as on matter degrees of freedom. After the 
holonomy replacement, this equation reads 



sin^(/(p)c)7b| 8^G 

9 r-?r \ 1 5— ^matter — U , (4: 

picking the sine as a real combination to be specific. (The precise form of this function replacing c^, restricted only to 
be bounded as a linear combination of exp(±«/(p)c) and to approach for small f{p)c in a low-curvature regime, can 
be motivated from full quantizations of curvature [ill Isol . [sij but is not determined uniquely. Further quantization 
ambiguities arise at this stage in addition to those in the form of f{p)-) 

When quantized, Eq. Q, provides the difference equation of loop quantum cosmology. For general matter fields 
the resulting behavior can be hard to analyze, but in some special cases the dynamics becomes exactly solvable. This 
is the case in particular for a free, massless scalar with Hamiltonian i?mattor = ^P~^^'^p'^ in terms of the momentum 
Pd, con j ugate to 0. Quantizing the system, subject to a particular factor ordering then results in harmonic dynamics 
[3 21 l33j . which leads to an evolution that is not subject to quantum back-reaction; in the sense that while states 
may spread and deform as they evolve, the changing shape does not affect the motion of expectation values. Hence 
quantum geometry corrections are dominant, resulting in the quantum bounce and the boundedness results. 

Harmonicity is realized thanks to a linear algebra of a complete set of operators including the Hamiltonian in 
deparameterized form. We thus first use the free scalar as an internal time variable, formulating evolution not with 
respect to proper time, so far appearing in derivatives denoted by the dot, but with respect to (j). This requires 
the monotonicity of (pit), which is globally guaranteed in some special cases including a constant (or vanishing) 
scalar potential, which are the cases considered below. Evolution for the relational quantities p((/)) and c((/)) is then 
generated by the Hamiltonian p^{c,p) obtained as a phase-space function by solving (|4]) for in -ffmattcr- We obtain 
P4, = ±^167rG/3(l - a;)|ImJ| with 

J:=Fexp(i/op"c). (5) 

The Hamiltonian is linear in J, and J forms a linear algebra with V by taking Poisson brackets: 

{V, J} = -iJ , {V, Jj^iJ , {J, J} = 2iV . (6) 

Classical equations of motion are linear and easy to solve, but even the quantum dynamics based on the Hamiltonian 
H := |ImJ| (dropping the constant ^/l6^^G/3{l — x) for simplicity, i.e. rescaling our "time" (j)) and the commutator 
algebra 

[V, J]^hj , [V, .P] = -hp , [J, P] = -2h{v + h/2) (7) 
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simplifies strongly.-'^ For instance, equations of motion d{V)/d(j) = {[V,H])/ih and d{j)/d(f) = {[j,H])/ih for ex- 
pectation values are fully determined by the expectation values and remain linear.^ This feature, shared with the 
harmonic oscillator in quantum mechanics, implies the absence of quantum back-reaction i.e. expectation values evolve 
irrespective of other moments of a state. 



B. Small potential 

Harmonic models are very simple but also extremely special. For robust conclusions, one must at least see how 
slight deviations from solubility might change the conclusions derived in this simple case. In the present context, the 
main ingredient for realistic modeling is the inclusion of a non-trivial matter potential. A number of key features 
change when the scalar potential no longer vanishes. First, cj) can at best be taken as a local internal time since it 
may develop turning points in a varying potential. Secondly, the Hamiltonian is no longer in a linear algebra with 
basic operators, equations of motion for expectation values no longer form a closed set, and quantum back-reaction 
ensues. Dynamics can no longer be formulated in terms of expectation values alone; we must include at least a certain 
number of moments of the evolving state, which, following |17|, we define as 

V---V J - J 

^ {{V {V)r{J - {J))')w.yi (8) 

for all positive integers a and b such that a + b > 2, indicating the totally symmetric ordering of operators by the 
subscript "Weyl". In general, all these moments are coupled dynamically to the expectation values, providing a 
complicated infinite-dimensional dynamical system. Approximately, it becomes tractable when an initial semiclassical 
state is involved, whose higher moments are suppressed by powers of h}^^ compared to lower ones. We will make use 
of this approximation in our analysis. 

In the presence of a non- vanishing potential W{(f>), the expression for on the constraint surface satisfying 

AttGpI vol^/^(Im e'-^^?')")^ 

— -\- vo\W = — 

3 vol 7^/(p)^ 



If we define 



3 

and absorb the same factor in as in the free case, we see that the momentum is proportional to 



H{V, J) := y^(ImJ)2-F3/(i-:.)p^(0) ^ (10) 

constituting a non-linear and typically "time" -dependent Hamiltonian. The rescaled potential W depends on Vq (via 
/o) in such a way that ^/^/(i-x)^ behaves like Vq. 

We will assume the potential to be always small compared to the kinetic energy of the scalar, putting us in the 
kinetic-dominated regime of loop quantum cosmology, close to the harmonic model. Sine the harmonic model has a 



^ Changing Vq is subtle at the quantum level. Classically, the Poisson bracket is not preserved, implying that the transformation is not 
canonical. The non-canonical classical transformation then is not implemented unitarily. One can mimic the change of the Poisson 
bracket by a rescaling of h, which does not amount to a physical change of its value but rather of the representation. Observables even 
at the quantum level do not depend on Vb when these properties are taken into account. 

^ With the absolute value in H = |ImJ|, the Hamiltonian is not strictly linear. However, the Hamiltonian is preserved in the case of a 
</>-independent potential, and so one only has to ensure that an initial state is supported solely on the positive (or negative) part of 
the spectrum of Im J and the absolute value can be dropped for all the evolution. For a state re quir ed only to be semiclassical at one 
time, the spectral condition can easily be arranged. Alternatively, effective-constraint techniques pj4436l | exist which allow one to deal 
directly with the constraint, avoiding deparameterization and taking square roots. 
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large class of states that remain semiclassical (or even coherent) for long times even across the bounce, the kinetic- 
dominated semiclassical regime of anharmonic models can be expected to allow long evolution as well, an expectation 
to be confirmed in this article. Larger potentials, for which states and bounce properties can change more drastically, 
will require significantly more care. 

For a small potential, we can expand the square root, but terms containing the potential remain non-linear in V 
and J. Adding the potential to the system thus provides an anharmonicity, and the exact solubility properties of the 
free model disappear. Not only the classical but also the quantum dynamics become much more complicated and 
in general can be treated only by using approximations. Even if it is possible to define the quantum dynamics in 
deparameterized form, such as with a constant W as used below, we have to face the difficulty of defining a complicated 
square-root operator for p^. 

At this stage, effective equations offer valuable tools to approximate the quantum evolution of general semiclassical 
states, sidestepping many difficulties. For a semiclassical state, we can expand the quantum Hamiltonian 

HgiV, J, G^'") := {HiV, J)) = {H{{V) + {V - {V)), (J) + ( J - (J)))) (11) 

as a Hamiltonian function on the quantum phase space of states parameterized by the expectation values V {V) 
and J :— (J) together with the moments G"'^. All the moments, as well as the expectation values, are subject to 
Poisson bracket relationships that are uniquely defined by {(A), (B)} = {[A,B]) /ih\ thus, the quantum Hamiltonian 
generates Hamiltonian equations of motion AF/A(f){y^ J, G"''') — {F, Hq} as usual. Also in this expansion the quantum 
Hamiltonian is clearly non-linear for a non-vanishing potential. Coupling terms between expectation values and 
moments imply quantum back-reaction. For semiclassical states, high moments are sub-dominant and the infinite 
sum can be truncated to finitely many terms, starting with quantum corrections due to fluctuations and correlations. 

We note that our condition for semiclassical states is rather weak, and so the allowed class is very general, much 
more general than usually considered in evolutions of specific states. We only require that a moment of order a -I- 6 -I- 1 
is suppressed by an additional power of h}^^ compared to moments of order a + b. (For dimensional reasons, the 
behavior must then be C''' ~ (VH)^'^'^^^^'^ in our variables, V and J having the same dimension as h.) This behavior 
is realized for Gaussian states (which fix all moments in terms of just one or two real parameters), but certainly not 
only in that case. States are allowed to spread out, deform their shape and deviate from Gaussians. As long as the 
hierarchy of moments is satisfied, the approximate quantum evolution remains reliable. In quantum cosmology, the 
evolution equations we use can hold true for sufficiently long times to answer long-term questions, and their validity 
can be tested self-consistently. 

In the specific case at hand, the quantum Hamiltonian is 



^ 2i J- J 2 (l-x 2 J- J 



J -J V^3/(i--) _ 3 2 + x l/(i+2x)/(i--) 

i 

h J- J 

3j y{2+x}/{1~x) 

l-x {J-Jf 

^3/(l-j;) 



(G^-^ -G^^)VK(0) 



(J- J)3 



[G-'-^ - 2G'^'^ + G'^'')W{(I)) (13) 



to second order in the moments, combined with an expansion up to linear terms in the potential. Note that the energy 
of a system of physical volume, vol, is related to Hq as 

1 SttG 

Energy = --^ -f vol = —(1 - xfHl/Yo\ + vol , (14) 

so Hq I vol is proportional to the square-root of the kinetic energy density. 

In the following, V and J refer to expectation values of the basic operators in a state. Their products in interaction 
terms of the Hamiltonian thus are not subject to factor ordering ambiguities as they would if Hq were an operator. 
Also in the underlying Hamiltonian operator there are no ordering ambiguities since the expansion comes from a 
Hamiltonian quantizing (jlOp . which lacks products of non-commuting operators. Although we have not written down 
an operator quantizing the square root in ()10p explicitly, it is clear that the semiclassical expansion is unaffected by 
ordering choices. This statement holds provided that we choose J as one of our basic operators together with V, as 
is done for the free harmonic model to make that system's Hamiltonian linear. Choosing J as the basic variable as 
(opposed to, say, f{p)c) is consistent with our use of moments of V and J. 
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We note that many quantum effects can be summarized as an effective Friedmann equation |37l . |38| 

(1)'.'^ L (: - i5L) , i, - P) + '-i^A (15) 

\aj 3 V V Petit J 2V Pcrit p + P J 

witli tlie critical energy density pcrit — S/SttGj'^L'^ and quantum corrections from correlations rj as well as fluctuations 
entering pg (and P is pressure). However, quantum correlations and fluctuations are determined by the moments of a 
state, and those moments are dynamical. An effective Friedmann equation thus provides a consistent set of differential 
equations only if it is accompanied by evolution equations for the moments. In this article, we are dealing from the 
outset with all equations for expectation values and moments to a certain order, bringing all quantum degrees of 
freedom under control. 



C. Non-dynamical conditions 



In addition to the dynamics of evolving moments, further conditions arise. First, to realize solubility of the potential- 
free system we have chosen to work with a complex variable J, which classically must satisfy | Jp = for the original 
phase-space variables c and p to be real. Reality conditions must also be imposed at the quantum level so that real 
phase-space variables are promoted to self-adjoint operators, ensuring that one is using the correct physical inner 
product to normalize states and compute expectation values and moments. After quantization, the reality condition 
still holds but is quantum corrected; we have 

|J|2_(T/+ lft)2^G^^_G^.'>in2 (^g) 

from taking an expectation value of the quantum reality condition 

Jp = (17) 

and expressing this in terms of moments and expectation values. The reality condition thus relates second-order 
moments to expectation values, but unlike the quantum Hamiltonian it is not truncated: it is an exact relationship 
for a finite number of our quantum variables. Note that in this expression one can neglect the terms, because we 
work in the regime V ^ h and the second order moments are order Vh and hence dominate the right hand side. 
However for clarity we will keep these terms and take the appropriate limits only at the end of calculations. 

In addition to (|16p . which relates expectation values in a semiclassical state with second-order moments of order 
h, reality conditions exist for the moments themselves. By taking the expectation values of Eq. (jl7|) . one arrives at 
Eq. (jl6p . however one can find further independent conditions, implied by taking an expectation value for instance of 
V^Jjt = . These new conditions relate third-order moments 

. ((f-(f))(i-(i))(jt-(Jt)))^^^^ (18) 

G^^^ ^ {{v-iV))') (19) 

to second-order ones and expectation values. Our truncation considers the evolution equations only up to second- 
order moments, disregarding third-order or higher ones, nevertheless, reality conditions obtained at third order restrict 
allowed choices for second-order moments, and thus will be important irrespective of the precise values of third-order 
moments. For instance, we obtain 

^VJJ _ qVVV ^ (^y ^ ^ qVV _ 2 (Re G^^Re J + Im G^'^Im j) - ^ [v'^ -^y , (20) 

an equation which contains terms of different orders in fi. Products such as VG^^ are of order in a semiclassical 
state, while third-order moments behave like {yKfl"^ and are thus smaller. The leading order in K then implies a 
reality condition 

VG^^ = Re( JG^^) + O (fr''/^) , 

= Re JRe G^'^ + Im Jim G^'^ + O ({VKf^'^^ , (21) 
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for second-order moments. By analogous calculations, we obtain 

VG^-' = ^{JG'''' + JG^-^) + O (^{Vhf/^^ (22) 

from the leading order of an expectation value of = JV'^. Eq. is a complex equation, giving two real 

conditions 

VRe G^-^ = i (rc JRe G" + Im Jim G^"' + Re JG''-^^ + O ([Yhf^^ , 

yim G^^ = i (rc Jim G" - Im JRe G'^"' + Im JG"'"^) + O {{Vhf/'^^ . (23) 

Similarly, low-order moments arise in expectation values of products with even more factors multiplying J ,V — V^.^ 
One could worry that arbitrarily high orders must be considered in order to find all reality conditions at second 
order. Fortunately, a combinatorial argument shows that this is not the case, and the reality conditions found here 
for second-order moments are complete. This conclusion can also be supported by a simple counting of degrees 
of freedom: We expect three independent second-order moments once all reality conditions are implemented; two 
fluctuations and one correlation parameter. Counting real and imaginary parts separately, we start with six moments 
of the (V, J)-variables (G^^ and G'^'' are always real, in contrast to G^'^ and G'^'^). With three reality conditions 
for the moments, a real one (j2ip and one complex one (1221) . the correct number results. 

In addition to reality conditions, moments of a state must satisfy uncertainty relations. In our case of a non- 
canonical set of basic operators, there are three independent conditions for the second-order moments, which can be 
derived by standard means: 

2G^^ (Re G''^ + G'''^^ - 4(Re G^'^)^ > h^ilm jf , (24) 

2G^^{-ReG''-' + G-'-')-A{luvG^-' f > h^iRe J)^ , (25) 
(G-^^y - (Re G-'-')^ - (Im G'^^)^ > h^iV + . (26) 

Note that for near-saturation these equations are consistent with the behavior G'^'^ ^ Vh of second-order moments 
with a + b — 2. 



D. Evolution of the Moments 

Following we obtain equations of motion for all the moments, in addition to quantum-corrected equations for the 
expectation values, by taking Poisson brackets with the quantum Hamiltonian. These equations, even when truncated 
to second order in the moments, are long, and for that reason given only in Appendix |X1 but as a finite-dimensional 
dynamical system they are amenable to analytical as well as numerical investigations. The most important one of the 
equations is the evolution of V, since this will determine the conditions under which a bounce occurs; it is given by 

6 J±J_y^2+.)/(i-.)(^^vj _ G^J)wi<p) 

[G^-^ + G^'^)W{(j)) 



l-x{J -J) 

3 y(2+2;)/(l-x) 

'l~X [J-Jf 
3 2+x J+J 
^2{l-x)^ {J-Jf 



y{l+2x)/{l-x)QVV ^^^^^^ 



-H3 ^ + ^ T/3/(i-x)(^.7j ^ gjj _ 2G-'-^)W{c^) . (27) 
(J — J) 



Interpreting J J — = as a constraint on the classical phase space, this large set of expectation values is an example of effective- 
constraint methods developed in [3343611 . The reality condition, although just one constraint, is of second class on the non-symplectic 
phase space spanned by {V, J, J), using generalizations of Dirac's classification of constraints to general Poisson manifolds 39]. Thus, 
only conditions for the moments result but no gauge flow need be factored out. 
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In total there are nine evolution equations for the real variables 

(v, Re J, Im J, G^^, G'"\ Re G^-^ , Im G^'\ Re G"^"', Im G^'') , (28) 

but not all the variables are independent owing to the non-dynamical conditions: the reality condition (jl6l) combined 
with restrictions from the uncertainty relations. Moreover, ([2T|) and ([22l) . while not presenting sharp conditions 
when third-order moments are not specified, provide ranges allowed for the moments of an initial semiclassical state. 
Imposing the reality condition (jl6p leaves eight free variables. In addition, we will require that the uncertainty 
conditions are nearly saturated at least for an initial state. As a sharp condition, this can be used to fix the initial 
values of, e.g., Re G^"^ , Im G^'^ and G'^'^, while the reality condition can be used to eliminate (for example) Im J, thus 
leaving five variables: V, Re J, G^^, Re G'^'^ and Im G'^'^. The relationship for the deparameterized Hamiltonian, 
Hq oc constitutes another equation which could be used to eliminate one further variable in favor of p^. (But note 
that is a constant of motion only in the case of a constant potential; otherwise it is subject to its own, potentially 
complicated dynamics.) The remaining moments then constitute two fiuctuation parameters, G^^ and Re G'^'-' , and 
one correlation parameter Im G'^'^.^ For numerical evolution, we will choose initial values for these variables and 
ensure that this is compatible with (PT|) and (f^ . 

Of particular interest is the evolution of the correlation parameter Im G'^'^ because its value plays a role in seeing 
whether the bounce of the free system generically remains if it is perturbed by the potential. Moreover, it contributes 
a term to the effective Friedmann equation ([T5|) similar to a positive effective potential [s^, HI] . Effective equations 
are not yet available in complete form for a non-trivial potential W{(f>), in whose presence </> would not serve as global 
internal time. Additional corrections can be expected in that case, but they can be shown to be small for a sufficiently 
fiat potential [40] . In our numerical discussions below (Sec.|V]), we avoid this issue altogether by working with constant 
(but non-zero) potentials, sufficiently small to ensure kinetic domination as used in the expansions. Our goals are 
to provide a more systematic analysis of consequences of potential terms compared to what is available so far. In 
particular, we will analyze the effects of a non-zero potential on the appearance of bounces and the evolution of the 
moments. 



III. ANALYTIC RESTRICTIONS 



Typically one is concerned with a state that is semi-classical initially, which in this language means that the moments 
are small compared to the expectation value. It is useful then to define the fractional moments, 

and analogously for higher-order moments AG"'''. Fractional moments in a semiclassical state then behave as AG"'** ~ 

(;i/V^)(a+6)/2. 

Note that the second-order moments AG^^ and AG"''' (but only those) must be positive. In the free model, V and 
I J| > ImJ cx take macroscopic values at the bounce for large matter content in the kinetic-dominated regime, all 
fractional moments should thus be small if an evolving state remains semiclassical well into the bounce regime. The 
state remains semiclassical at all times in the free model, so that a sufficiently small anharmonicity brought about by 
the non-vanishing potential can be expected to allow semiclassical evolution into the bounce regime for a reasonably 
large class of states. 

There are additional restrictions on specific combinations of the moments, imposed by uncertainty relations. In 
terms of the fractional moments, those relations, Eqs. (f24|) -([26 l) . become, 



2AG^^ (Re AG-^'^ + AG^'^") - 4 (Re AG^'^' > ^^^^ ' (^0) 
2AG^^ f-Re AG'" + AG'"1 - 4 (im AG^'^)' > h^''^^^ , (31) 



* To classify second-order quantum variables as fluctuations or correlations, it is preferable to refer to the self-adjoint operators J+ : = 
^(J -I- J^) and J— := -^{J — J^) rather than J and jt. Since this is a linear transformation, the order of moments is preserved. In 
particular, G'^'^ = G'^+ — G''^ ''^ + 2jG'^+ ''^ , wrhose real part has only products of the same operator, thus constituting a fluctuation, 
while the imaginary part contains the product of two different operators J± . 
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(ag'-^Y - (Re AG^'Y - (Im AG^'f > fi^IKl^^ ^ (32) 
V / \J\ 

Using Eq. (|16p and Eq. p2p we find an inequality for V, 

h 



V > 



22?(1 + AG^^) 



1 + AG'''' - V 



l + AG-J'^ - vf + 2A (1 + AG^^) (1 + AG'^-^" - V) 



where 



(33) 



V=\l {AG-^Jy - (Re AG^-^)' - (Im AG^^)^ . (34) 
In addition, from Eq. ([5(11) and Eq. pip , we find, 

V > ^ . (35) 

2\l AG^^AG-" - (ARe G^'f - (Aim G^'f 

These two conditions must always be satisfied, and in particular they must be valid at the bounce, which is by 
definition the minimum value of the volume. Via Eq. ^ we find that the volume at the bounce is restricted. In 
particular, Eq. psp implies that, 

(aG^^ AG^'^" - (Re AG^'^) ' - (im AG^'^) ^ , (36) 

V / Bounce 

can only occur for Vbounco — > oo. Thus a semi-classical state (at the bounce) is not just one in which the fractional 
moments are small, since the uncertainty relations for such a state, would force the volume (of any spatial region at 
the time of the bounce) to be arbitrarily large. Second-order moments thus cannot be self-consistently ignored. 



A. Zero potential 



In this section we consider, in more detail, the case of a free scalar field, i.e. W = leaving the discussion of the 
W ^ case to Section ITlI B I The evolution of V is given by Eq. (P7)) and in particular the bounce is given by ^ = 
which gives. 



Re J 



. 



Putting this into Eq. (jl6p we find. 



Im J 



V^l + AG^^)+n{V + ^) 



-I 1/2 



(1 + AG") 

where the second equality comes from Eq. (jl3p . Thus we find the density. 



Hq 



Hq 
vol 



vol 



.1 V {l + AGVV) + n{V+\) 
V d + AG'") 



(37) 



(38) 



(39) 



It is important to note that although this resembles an energy density, we have Hq oc p^, and hence the energy 
density is the square of Eq. (15^1) (see the comments below Eq. P^ ). 

If we assume that at the bounce, V is large^, i.e. that ^ » ?i, then the density at the bounce is given by the 



^ The variable V depends on Vb and can be made arbitrarily small (provided a; < 1) just by choosing a small Vq. But we have to be more 
careful at the quantum level due to the non- unitary implementation of changing Vq. For simplicity, we consider the limit of large V as 
a substitute of the complete transformation, and ensure that physical values remain finite in the limit. (The limit of large Vb has also 
been discussed in [4ll .) 
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large- limit, which is just, 



Hq 

vol 



Bounce, larec V 



SnG-f (l-x) 



(40) 



where the fact that (l + AG^^) > 1 and yl + AG'-''-' j > 1, ensures that the square root is real. Note that, as 

expected, this is independent of the fiducial volume Vq (and f^^p^(^+^^y^ — L^^ is constant in our limit of large V^) 
and so is the bounce density 

AG^^ 1 + AG^^ 



P 



_ 3 1 

Bounce, large V 



Pcrit 



1 + AG-^-'" 



(41) 



The p- independence of the bounce density for the specific case of a: = —1/2 agrees with numerical results [42|, as 
pointed out in [s^] whose more general results are consistent also with the other cases of x. The specific dependence 
of the bounce density on the moments agrees with results obtained in the context of the effective Friedmann equation 
(ITSI) . where the energy density pq = p + copait + ■ • •, according to [13, HI], has the leading correction given by 
eq = AG^'-' — AG^^ . There is a bounce when pq — pcrit, at which the density is p = Pcrit(l — eo)- For small fractional 
moments, this is exactly (PT|). 

Notice that even for the x — —^^ case, the energy density at the bounce still depends on the moments present at 
the bounce (even with = 0, i.e. zero potential). The presence of these moments can increase or decrease the energy 
density of the bounce compared to the case when they are neglected, however for semi-classical states, we can restrict 
the energy density to the range, 



Bounce, large V 



< V2p. 



'crit 



(42) 



where pcrit is the energy density of the bounce that would be calculated neglecting the effects of the moments and we 
consider a state to be semi-classical only if the factional moments are less than unity (admittedly, a rough estimate) . 
The cases of equality can be removed by noting that Eq. (15(11) implies AG^^ > and Eq. (15^ gives AG''"' > 0. This 
is a limit on all possible semi-classical states, for which the higher order moments can be neglected. If one wanted to 
consider a more strict definition of semi-classical states (i.e. states for which the fractional moments are below some 
specified value, less than one), the range would be further restricted. 



B. Non-zero potential 

For the more interesting case oi W ^ 0, Eq. (P7| and the bounce condition dV/d(t) — gives, 

WV^ f 1 Re J 3 Re J IJP 



Re J = -- 



Im J 1 2 Im J 4 Im J (Im jy 



Re AG'''^ - AG 



3 Re J I J| 
1 - a; Im J (Im J) 

1 IJI2 



rim AG 



VJ 



rim AG 



7.7 



3 2 + x Re J 
4(l-a;)^Ini J" 

1^1 



AG 



vv 



2{lmjf"" ' 2(l-x)ImJ 

Changing to polar coordinates for the complex variable J using. 



Re AG 



VJ 



Im J 



smf = 



where 9 — arg (J), we can write Eq. (|43l) as 



cos 9 



WV 



cost 



1-== r cost 
in0l2si 



Re J 



3cos( 



tan 9 — 



Im J 
Re J ' 



J|2 sin 6* I 2 sin 6* 4sin^ 



[Re AG^^ - AGjj] 



ImAG^^-H H2 + x) cos9 ^^yy 



(1 — x) sin^ 
1 



2 sin^ 9 



Im AG 



4 (1 — x) sin t 
^ Re AG 



VJ 



2 (l-x) sin 6* 



(43) 



(44) 



(45) 
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Now we look for solutions that are 'close' to the W = case in the near-bounce regime, by expanding around 
9 = tt/2 i.e. = n/2 + 60, which gives sin6' = 1 + O (69^) and cos 6* = -69 + O (69^). We will see that this 
approximation is valid within the accuracy of the expansion. With Eq. (fT6|) we get, 



^Re AG^-^ - ilm AG-'^ 



Again, taking the large-V limit, we find 



^jT^f^Re AG^-^ - ilm AG^'' 

2(1 — x) 2 



large V 



1 

2 4 



^ ^Re AG-^^ - AG^-n - v^Im AG^^ + fy^i^AG^^ 

/ L—X 4 [1 — x) 



(46) 



Thus the expansion is valid for small fractional moments: 69 is small if Aim G'^"^ and ARe G^"' are small. In this 
approximation, there is always a solution to the bounce condition (j43p . It is clear that the bounce cannot disappear 
as long as this approximation is valid. We are thus not testing the robustness of the bounce itself, but rather some of 
its properties, such as its energy density. 

We can now look at the (square-root of the kinetic) energy at the bounce which from (|13p is given by. 



Hr 



\J\ sin I 



W 



Bounce 



|J|sin6' 



1 , 3i'^ + ^)_^^vv _ 3ImAG^'^ 



1^2 4(l-a;)' 
Once again, expanding about 9 = Tr/2 + 60, gives 



2(1 -a;) sine* 4sin2 6' 



Re AG"'-'' - AG 



-.J.7 



(47) 



Hr 



\J\ 



1 3(2-|-a;) 



4(1-2;)' 



AG^^- 



VJ 



3Im AG 

2(1--^) 



4( 



Re AG 



jj 



AG 



JJ 



o{60y 



(48) 



Here it is important to remember our redefinition of the potential: every W {(j>) comes with a Vb-dependent factor via 
/o in ([9]). Using this and the reality condition Eq. (fT6|) one finds that the (square-root of kinetic) energy density at 
the bounce, again taking the large-F limit, is given by. 



Hq 
vol 



Bounce, large V 



A~WA 



1 + AG" \ ( 87r7G/o2pi+2^ ^ 



1 + AG^^ j \ 



1 . 3(2 + ..)^^^,,_3ImAG^^_l. 



4(1-2;)^ 



where 



A 



2(1 -x) 



mi 

vol J large y 



O{60f 



(49) 



(50) 



Using the reality condition Eq. ([TB)) and Eq. ^ we see that this is 

3 



A = 



BttGi (1 - x) 



1 + AG^v 
1 + AG-^-^" 



r(0) , 



(51) 



which is just the (square- root of the) energy density at the bounce that we had for the W = case, Eq. (15^1) . So 
finally we have that while the bounce is 'close' to the W = case, its (square-root of the kinetic) energy density is 
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given by, 



r{W) = r{0)-Wr{0) 



3{2 + x) A^^yy 
4(l-x)' 



— -Im AG^^ + l(Re AG^-^ - AG"' 
2 (1 — a;) 4 



o [sey 



where 



Bounce, large V 



(52) 



(53) 



Just as in the W = case this is exphcitly independent of the fiducial volume. 
For the x = —1/2 case, for instance, this is 



1 1 



-AG 



vv 



Im AG' 



Re AG 



jj 



AG 



J J 



-o{5ey 



and the small parameter that we have expanded in is, 

ARe G^^ 

large V 



(54) 



59 



ilm AG'^"' 



1 _ 3 

2 4 



(ARe G'J-J - AG'^'^") - 2Im AG^'^ + iAG^^ 



In general, we can say that the (square- root of the kinetic) energy density at the bounce for x = —1/2, for semi- 
classical states (i.e. states with fractional moments less than unity), is restricted to the range. 



V2r(0) 



^7rG7/o^ 



I 1 + AG^^ j 



>r{W)> 



V2 



1 + AG 



with 



/ = l + ^AG^^ 
2 2 



Im AG^'^ - - f Re AG'''^ - AG'^-^" 
4 



(55) 



(56) 



where we used Eq. and, as before, r(0) is the (square-root of the kinetic) energy density that would have been 
calculated neglecting the moments entirely. 

To decide whether the effect of moments raises or lowers bounce densities, it is more convenient to consider the 
total rather than kinetic energy density. After all, from the effective Friedmann equation ([15]) we expect the total 
density to trigger bounces. With (jl4p . we have the energy density at the bounce given by 



SttG , \9 / Hq 
PBounco = — j^(l - a;) 1-^ 



1 -<r AG^^ 
W ^ pcrit . -DW 



1 + AG-'"' 

with the critical density pcrit expected for the free bounce, and (for x = —1/2) 

1 



D = 2/ - 1 = AG 



vv 



2Im AG 



VJ 



Re AG 



J J 



AG 



J J 



(57) 



(58) 



Also this correction to the bounce density agrees with the term 5i from [33| (or ti from (Sgl) based on the effective 
Friedmann equation. By Eq. (j31[) we have 



(Re AG"'-' - AG'-') > 



(59) 



and AG^^ is positive while Aim G^"' does not have definite sign. The sign of D can thus take both values, raising 
or lowering the bounce density on top of the effect due to the potential.^ 



In [43l . it was shown that expectation values of a density operator, of the form (Energy/volume), are bounded from below by the 
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IV. DYNAMICAL COHERENT STATES 



Of particular interest are states which saturate the uncertainty bounds Eqs. p4l)- ((26l) and remain saturating them 
throughout the evolution of the system. Such states are referred to as dynamical coherent states. We will first 
construct such states for the W {(/)) = case which is analytically tractable and then discuss the consequences for the 
rather more complicated case of W {(/)) ^ 0. 

For a state to initially saturate the uncertainty conditions we require, 

2G^^ (rb G'^'^ + G'^'^") - 4(Re G^'^f - ^^(Im J)^ = , (60) 

2G^^(-Re G'^'^ + G'^'^) - 4(Im G^"' f - ^^(Re J)' = , (61) 
{G-^^f - (Re G^^)2 - (Im G^-')'^ -h^{V + h/2f = . (62) 

These states must also satisfy the reality condition, Eq. (IT51) . which we rewrite here for convenience, 

Dynamical coherent states (as with any other state), must also satisfy Eqs. (|2i p -(|22 |) . however these are 'weak' 
constraints in situations in which third-order moments are not specified, in the sense that they need only be satisfied 
up to order iyhY/"^ and hence do not provide sharp constraints on the states. 

A state that satisfies these four conditions will initially saturate the uncertainty bounds, however to ensure that 
it remains saturating them we also require that the derivatives of Eqs. ([60 |) -(j62 |) with respect to be zero. For the 
W ((/)) = case one can easily show that the derivative of Eq. (1501) with respect to automatically vanishes, while for 
Eq. ((6T|) and Eq. ([62)) we require, 



2Re G 



VJ 



Re G" - G" 



2Im G^-'lm G^-^ + ft^ ( F + ^ ) Re J = . (64) 



This one constraint is sufficient to ensure that both Eq. (|6ip and Eq. ([5^ are independent of cf). This makes the space 
of dynamical coherent states a hypersurface of co-dimension one in the space of states that instantaneously saturate 
the uncertainty conditions. Explicit solutions have been found and discussed in [l^,!!!]. The existence of a large class 
of coherent states is analogous to the well-known feature of the harmonic oscillator in quantum mechanics, but notice 
that potentially significant spreading of states is generically possible even for dynamical coherent states [1^ [20] . 

The full parameter space of this system is 9-dimensional, however the five constraints above can be solved to give 
dynamical coherent states in a 4-dimensional parameter sub-space. A particularly convenient way of solving these 
constraints is to specify initial data 

ft/, Im J, Re G", G'^'^) , (65) 



solve Eqs. ([60])-([62|) as 



Re G^^ = ±\/ \g^^ (Re G-^J + G^J) - hi^{lTn Jf , (66) 



Im G^^ = ±\/iG^^(-Re G-" + G^-^ - ^^^^(Re J)2 = , (67) 



Im G'' = ±V(G-^-^)2 - (Re G^-^f - h^V + h/2f , (68) 



critical density of the harmonic model even in the presence of a cosmological constant. Our statement is not in confiict with this 
result because we are considering bounce values for expressions of the form (Energy) /(volume). Obviously, the two expressions differ 
by terms depending on the state via its moments, exactly the kind of terms we compute; what we show is that those moments may 
pull the bounce density below the critical one of the harmonic model. In this context, one should note that l43l refers to a density 
operator which cannot exist in full loop quantum gravity where only total Hamiltonians, but not energy densities of matter or the 
gravitational field can be represented as well-defined operators. Expressions such as (Energy) / (volume) used here, on the other hand, 
have completely well-defined analogs in the full theory. Thus, our results have a higher degree of robustness than those of |43|| . even 
though the mathematical statements of [43l | are sharper. 
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and then simultaneously solve Eq. (163]) and Eq. (|64)) for and Re J. There is an analytic solution to this 

system, however in practice it is more useful to solve the system for some given initial data i.e. to specify 

l^y, Im J, Re G'" ^ G'^"^^ and solve the resulting simultaneous equations. This is a useful test of the numerical 
implementation of the system, which we discuss in the following section. 

In order to attempt a construction of dynamical coherent states for the W {(/)) 7^ case, one can in principle follow 
the same procedure. Clearly the initial states will still have to satisfy the constraints given by Eqs. (|60p - (p5)) . however 
in order for these conditions to be independent of (j) we now need to additionally satisfy. 



,dG^ 
' dcj) 



Re G'-' + G^'^ 



2G 



vv 



2^. ReG^^ + G-^-' +2G^^ 



dRe G-'-^ 


dG" 


d(p ^ 
dRe G"'"' 


dcp 
dG-'^ 


d(t> ^ 
2G-^^ 


d(j) 
dG"' 



8Re G-^^^ 2nHm = , (69) 



-8ImG-^i^-2.^ReJ^ 



2Im G 



J J 



dim G 



JJ 



2n^ {v + 



h\ dv 



, (70) 



(71) 



where the derivatives of G""'^ are given in Appendix [A] These three equations need not be degenerate and certainly 
are not automatically zero for general states and general W {(j>). They correspond to the single equation Eq. (IM)) in 
the W {(j)) = case. These equations are significantly more difficult to solve. Naively it would appear that we now 
have seven constraints on the nine-dimensional parameter space that give dynamical coherent states for W (0) ^ 0. 
In fact it is easy to see that Eqs. (|M | - ([7T|) contain at least four independent conditions on the moments. To see this 
note that all of the derivatives given in Appendix [X] are of the form. 



TV — 



W (0) Fi 



(72) 



where Fo and Fi have no dependence on W {4>). It is then clear that Eqs. ([Mjl - dTTt contain terms that are independent 
of W {(p) and those that are not. For a general potential, these must independently cancel which gives us four 
conditions. The W {4>) independent part of Eqs. ([55)) - ([7T|) gives exactly Eq. (|64p . while the W {cj)) dependent parts 
give three additional constraints. Thus we are left with (at most) a one dimensional parameter space that can support 
dynamical coherent states for a general W {4>) . 

Thus far we have been concerned only with the conditions imposed on the moments by the requirement that the 
uncertainty conditions be saturated and their first derivatives be zero. This alone is not sufficient to fix a dynamical 
coherent state, since we have to require that all derivatives of the uncertainty conditions vanish. For the W {(j>) = 
case, one can explicitly check that the derivative, with respect to 0, of Eq. (|64p vanishes provided the uncertainty 
conditions are saturated (in particular provided Eq. (|5T|) and Eq. (|62p are saturated, since Eq. ([SO)) decouples from 
the system in this case). Thus by inspection Eq. ([M|) is the only non-trivial condition that needs to be satisfied. 

For the W {(/)) 7^ case we have four constraints to investigate. The W ((/>)-independent part of Eqs. ([551) - ([7T|) is 
just Eq. ((64|) and hence the previous argument holds. For the W {(p) dependent parts of Eqs. ((69)) - ([7T|) the situation is 
less clear, due to the complexity of the evolution equations given in Appendix[Xl However if any one of the derivatives 
of any of W {4>) dependent parts of Eqs. ([551) ~ ([7T|) does not vanish then the parameter space that supports dynamical 
coherent states reduces to a single point and if any two of the derivatives fail to vanish dynamical coherent states 
cannot exist. 

In addition to the fact that dynamical coherent state may not exist, this reduction in the dimension of the parameter 
space for such states is also a difficulty for study of the dynamics. Since one can no longer explicitly solve the eight 
constraints implied by Eqs. (pn|) - (|55)) and Eqs. (|M)) - ([7T|) ^. small changes in the initial state can lead to radically 
different trajectories. This is an inherent property of the evolution equations of the states and can thus play a 
significant role in both analytic and numerical investigations of the dynamics of (nearly) dynamical coherent states (if 
such states exist). Numerically, difficulties can arise from the large parameter space to be explored when one can no 
longer stay on the saturation surface. One may reduce the number of initial conditions by starting on the saturation 
surface, but numerical errors easily lead one into danger of falling onto the wrong side of saturation. 



Recall Eqs. I|69| l- H71|l are four independent equations. 
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V. NUMERICAL EVOLUTION OF THE SYSTEM 

The evolution of the moments and their back-reaction on the trajectories of (J, V) can be directly implemented 
numerically. In particular, because Eq. (I27p and the equations in Appendix |X] are flows in configuration space (i.e. 
depend explicitly only on the variables and not their derivatives) the numerical implementation is rather straightfor- 
ward. Care is however needed to ensure that our numerics are sufficiently accurate so as not to introduce violations 
of the uncertainty relations, Eqs. or reality constraints, Eq. (ITBl) . (PT|) and ([2^ . when exact conservation is 

expected analytically. In the W case, as already discussed, exact conservation is not expected since the absence of 
higher order moments are expected to introduce a fundamental limit to the accuracy of our evolution. Nevertheless, 
as long as the fractional moments are small, all conditions should be respected. 

Another possible point of concern for the general W {(p) ^ case, is that it is no longer guaranteed that (p will 
be monotonic. If it is not, then the deparametrization of the system (using (f> as an internal time) will no longer be 
valid. Such a break down would correspond to a non-unitary evolution of the system and is clearly not physical. This 
issue is easily avoided, by choosing the </> coordinate only in regions where it is monotonic and changing coordinates 
at turn-around points. In our numerical investigation, we avoid this subtlety entirely by focusing our attention on 
the W (0) = const ^ cases, for which </) is monotonic throughout the evolution. Here, we choose < so as to get 
into a bounce regime at larger densities and smaller volume than in the free case. 

In practice, fractional moments grow, which can result in the non-conservation of the uncertainty conditions and 
limiting the length of for which the numerical equations can be evolved. However this is sufficient to examine the 
bounce, for suitable chosen states. To monitor the growth of higher-order moments which are not included in the 
evolution equations, we will use plots of the higher-order terms in reality conditions: 

VG^^ ~Re{JG^^) =: V^Ri, (73) 

VRe G^-' - i (Rc JRe G" + Im Jim G'^'^ + Re JG^'^") V^=^i?2 , (74) 

VlwL G^"' - i (rc Jim G" - Im JRe G"^'^ + Im JG'''^) =: V^Ry, . (75) 

^From the general behavior of the moments, we expect R2 and R^ to be significantly smaller than the fractional 
second-order moments, e.g. J^i ^ (AG^^)'^/^ ^ {h/VY^'^ . In addition, it is rather difficult to find initial conditions 
that simultaneously satisfy all of the uncertainty conditions. Thus it is often more convenient to use the uncertainty 
conditions to fix some of the initial data, by requiring that the state under saturate Eqs. p4 p - (j26p by some constant 
i.e. that the state satisfy, for example, 

2 AG^^ (arc G'' + AG'"") - 4 ( ARe G^-^) ' = + ' ^^^^ 

and similarly for the remaining two uncertainty conditions. In general C/i, U2 and C/3 depend on 0, however, from 
Section HVl we see that for the W = Q case, Ui should be a constant, while U2 and should evolve similarly. 

For the numerical calculations we used units in which h = 1 and V is dimensionless. A typical example for W {(jy) — 
is given in Figure [TJ which had the following initial conditions: 

/ V Re J Im J AG^^ AG'" ARe G^-' Im AG^"' Re AG"'"' Im AG"'"' \ 
1^ 35.0 24.0 25.0 0.0331 0.0833 0.0377 0.0294 0.0250 0.0737 j ' 

with X = —1/2. These initial conditions under-saturate the uncertainty conditions, Eq. and also saturates 

their derivatives with respect to (Eq. (|M|)). In addition, one can explicitly check that Eqs. (I?T|) -(P ^ are satisfied, 
up to the correct order, throughout the evolution of the system. In the numerical evolutions, we have taken h — \, 
hence one needs to check that Eqs. (|211) - ((22|) are satisfied up to order -\/AGyy, (or other second order moments), 
since this is the order of the terms neglected. 

For this state the bounce occurs at V (^Bounce) = 25.6582, with the square-root of the dimensionless kinetic energy 
density r (0) = 0.974348. This can be compared to the expected value, given by Eq. (15^ . of r (0) = 0.974349, showing 
the excellent agreement of the numerical implementation with the analytic solutions. In the W (0) — case, the 
evolution equations given in Appendix |^ are exact and hence we can analytically calculate the (square-root of the) 
energy density at the bounce to arbitrary accuracy (for any given value of AGyy and AGjj at the bounce). This 
is not true for W {(p) ^ 0, since the evolution equations are valid only up to order O {W^). However even for the 
W = Q case numerical approximation of the derivatives introduces an error, this can become significant for states that 
(almost) exactly saturate the uncertainty bounds, where numerical artifacts may lead to unphysical trajectories. 

If the saturation of the uncertainty conditions, Eqs. (I5(H) - (|5^ . and their derivatives (Eq. (|64l) for the W {(p) = 
case) were met exactly then the state would saturate the uncertainty conditions and remain saturating them, through 
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FIG. 1: The trajectory V {(j)), evolved numerically for the initial conditions given by Eq. (|77p . for the zero potential {W (</>) — 0) 
case. Note that only the region around the bounce is plotted. This bounce occurs a,tV~ 25.6582 with a corresponding (square- 
root of the) energy density of 0.9743. Note that for this trajectory, the bounce occurs at (/>Bounce ^ 1.821815. 
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FIG. 2: The system was initialized with a state that both under-saturates the uncertainty conditions, Eq (|24p - (|26|l and 
ensures that the derivatives with respect to </> are all positive. Plotted are the values of Ui, U2, 4.0 x U3 for the W {(j)) = 
and W {(j)) — —0.05 cases. Recall that Uj, U2 and U3 positive indicates that the uncertainty conditions are under saturated 
i.e. physically acceptable. Clearly all the uncertainty conditions are satisfied, in particular at the bounce, which occurs at 
(j) ~ 1.8. Notice that for W = 0, Ui is essentially constant as expected from the discussion above Eq. (|64|) . Note that the initial 
conditions are deliberately taken so as to ensure that the system is far from the uncertainty conditions throughout the bounce. 



out the evolution in (j). For the state here, the initial data was chosen so that the uncertainty conditions, Eq. (j76p . 
were initially under-saturated and the the derivatives of these equations with respect to (f> are positive for both the 
W {(j)) — and W {(/)) = —0.05 cases. This ensures that, at least initially, the system is evolving towards the under- 
saturated side of the uncertainty condition hypersurface. Figure [2] shows how the uncertainty constraints evolve for 
both the = and W = -0.05 cases. 

For the case of W {(/)) = —0.05, the trajectory V {(p) is plotted in Fig. |31 with the initial conditions again given by 
Eq. ((77)) . In this case the kinetic energy density is found to be higher than the W {(p) =0 case, see Table U (which 
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FIG. 3: The trajectory V {(p), evolved numerically for the initial conditions given by Eq. (|77p . Again only the bounce region 
is plotted. The solid curve correspond to the W {(j)) — case and the dashed curve to the W {(j)) — —0.05 case. The bounce 
for this latter case occurs a,t V ~ 25.11053 with a corresponding square-root energy density of 1.00515. This is, approximately 
3.2% higher than the W {(j)) = case, a change mainly (but not solely) attributed to the negative potential energy; see below. 
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FIG. 4: The trajectories of all the variables over the bounce range, for the initial conditions given in Eq. (|77[). The solid lines 
are the W (</>) = case and the dashed lines correspond to the W (0) — —0.05 potential. From top left to bottom right the 
trajectories are: V (0), Re J (0), Im J (<;/>), AG^^ {cf>), AG''-^ (</>), ARe G^-' (<^), Im AG^-^ {(j)), ARe G''-' {(j)) and Im AG'^-' ((j)). 
Notice in particular that for W {<j)) = we have Im J a constant, which is expected since for this case we have Im J = Hq 
which is conserved and that several of the variables (Re J {(j)), AG'''' {(f>) etc.) evolve asymmetrically around the bounce point. 



Wi<l>) 





-0.05 


vol 


Rnnnrr 


0.97438 


1.01515 



TABLE I: The kinetic energy density of the bounce depends on the magnitude and form of the potential (W (0)) that is present. 
The initial conditions are given by Eq. (|77[) 
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Estimate 


Eq. ^,56'^ ^ 0.0005 


numerical 


vol 


Rnnnrn 


1.0038 


1.0052 



TABLE II: The energy density of the bounce for the W {(j)) = —0.05 case, calculated using the approximate analytic results 
Eq. (|48|) and Eq. (|49|) compared to the numerical result. 



gives the square-root of the energy densities). Using the values of the moments at the bounce, we can evaluate the 
expected Hamiltonian density from Eqs. and compare them to the numerical result, see Table HT] Notice 

that for this case the evolution equations given in Appendix [X] are only valid up to O (W^) = 2.5 x 10~^ and an 
accuracy beyond ~ 3 significant figures would require the use of higher orders in the expansion. We can see that 
the approximation made in order to solve Eq. (j54[) is accurate to, approximately, 0.5%. It is also worth noting that 
the moments, as seen in Figure 21 do not evolve symmetrically around the bounce point. In particular Re J, AG'^'^, 
ARe G^"' , Im AG^"' and Im AG"^"^ all decrease monotonically across the bounce (for this system). This may be 
significant for deparametrizing the system in order to use one of these variables as an internal time parameter, as 
suggested in [i^. 

For the case of constant W as studied here, the effect of varying the moments is small but observable (for the 
evolutions considered here). In Figure [5] the trajectories V {(p) for various initial AG"'"' in the range, AG'^'^jj^^^^^^ = 
(0.108, 0.189). For W {(f) = —0.05 the effect of varying the initial moments is even more pronounced than for W = Q. 
Figure [S] shows the dependencies of the (square-root of the) energy densities at the bounce on the initial value of 
AG^^ , for both the W {(f) = and W {(j)) = —0.05 cases. In the cases examined here, AG"^^ is approximately 
proportional to AG^^ at the bounce. Hence, by Eq. (^0]) . we would not expect there to be a significant dependence 
on AG''^ for the W = Q case, which is what is found. For the W ^ cases however there is a dependence, due to the 
terms proportional to W in Eq. (|48)) . 

In Section HVl we described how to explicitly construct dynamical coherent states for the W {(p) = 0.0 case. Given 
the initial data, 

(v, Im J, G''^, Re G"'"') = (105.0, 100.0, 125.0, -1.0) , (78) 

we can solve Eq. ([7S| . with Ui = U2 = U3 = 0.0 (which is required if the state is to be dynamically coherent) and 
Eq. ([M)) to find (for example) 

(Gyy, Re J, Re Gyj, Im Gyj, Im Gjj) w (82.52970, 32.985295, 51.15507, 70.195184, 67.035438) . (79) 

We have checked carefully that for this case (and others that are expected to be dynamically coherent) the uncertainty 
relations are indeed conserved up to numerical rounding errors (which can be made negligibly small by employing a 
suitably accurate discretization scheme). This provides a strong test of the accuracy of our numerical implementation 
Finally as an example of moments which evolve to break the uncertainty relations (in this particular case Eq. (161^ ) 
showing the breakdown of the approximation, in Fig. [7] we plot the evolution of Ui, U2 and for the initial conditions 

V Re J Im J AG^^ AG''-' ARe G^"^ Im AG^-'' Re AG"'-'' Im AG"'" 
35.0 24.0 25.0 0.0331 0.09159 0.04267 0.02363 0.0250 0.0737 

This initial state initially undersaturates all the uncertainty conditions {Ui ~ 18.1, U2 ~ 2623.9 and C/3 « 2103.6 at 
(j) = 1), but dynamically evolves to break the uncertainty bounds. As one can see from Fig. [7|the moments are no 
longer small for this state (AG"'"' and Im AG"'"' are larger than 10%), which indicates that the contributions to the 
dynamics from higher order moments can no longer be neglected. Note that even for this case, there is a bounce at 
(j) « 1.82, with square root of the energy density being pbounce ~ 0.97435 and despite the fact that the evolution of this 
state has significant contribution from higher order moments, this value is within 0.72% of that predicted by Eq. (|39p . 
Thus we see that our estimates for the (square root of the) energy density at the bounce are well approximated 
using only the first order moment, at least while the dynamics of the state do not breach the uncertainty conditions 
Eq. Once again we emphasize that this breaching of the uncertainty conditions is not an artifact of the 

numerical implementation, rather it is a genuine feature of the dynamics, due to the truncation at second order in 
the moments. 
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FIG. 5: The trajectory V {(p), evolved numerically for the initial conditions given by Eq. (|77p . Again only the bounce region 
is plotted. In this case, the potential was set to be — 0.05. The various plots correspond to different initial values of jj 

in the range AG-'-' = 0.107779 0.18926 (Gjj|initiai = 100.0 140.0, in steps of 5Gjj = 1.0). The sensitivity to 

bounce I 1 I , 

variations of the moments is significantly increased after the bounce, an observation extending the conclusions of [l9l. |20|]. 



VI. CONCLUSIONS 



Numerical investigations in loop quantum cosmology have so far only considered semi-classical states of Gaussian 
type, peaked around some particular value [l^jH^li^. The expectation value of such states has been shown to closely 
follow the classical trajectories, except at very high densities where quantum gravity effects, but primarily those of 
quantum geometry, produce a bounce. Implicit in these analyzes is the idea that a general semi-classical state is 
well characterized by its expectation values and small spreads alone and that the evolution of the spread and other 
moments does not affect the evolution of the peak. Here we have considered a moment expansion of a state around 
such exactly solvable Gaussian states as a way of probing the consequences of relaxing this assumption. In principle 
one should consider all orders of the moment expansion, however in order to make the system tractable, we have 
truncated the series at second order. Already this expands the configuration space from 4-dimensional {{V, 4>) and 
their momenta) to 18-dimensional (given by the eight independent variables in Eq. (j28l) . (j) and their momenta). 

The inclusion of moments introduces sets of non-trivial uncertainty relations that must be satisfied by the state 
at all points along its trajectory. Dynamical coherent states are those for which all the uncertainty relations are 
exactly saturated, throughout the state's evolution. We have shown that the conditions necessary to ensure a state is 
dynamical coherent, effectively reduces the dimension of the state's phase space i.e. the conditions form a constraint 
surface on which the states, and their trajectories, lie. For the zero potential case, there is only one condition that 
needs to be met to ensure that a state which initially saturates the uncertainty bounds will be a dynamical coherent 
state (i.e. will continue to saturate the bounds). For a system with a non-zero potential one has four independent 
constraints (in addition to the uncertainty bounds) which need to be met, just to ensure that the first derivatives of 
the uncertainty relations vanish. While requiring that the higher derivatives also vanish is likely to introduce even 
more constraints. This makes it extremely difficult (and probably impossible) to find exact dynamical coherent states 
for systems with a non-zero scalar field potential. While conceptually this may be a problem, especially for considering 
evolution from the infinite past (0 — >■ —oo) or cyclic models, in practice we are used to dealing with quantum states 
that are only approximately dynamically coherent, such as particle wave-packets in standard quantum mechanics. 
Provided the states are approximately dynamically coherent for a sufficiently long period of time (0) they can, for all 
practical purposes, be considered semi-classical. 

We have investigated the consequences of the moments on the presence of a bounce, for states that remain semi- 
classical (in the sense that all their moments are small compared to the expectation values and that the states obey 
the uncertainty conditions). In all cases the bounce was found to be present and the energy density was changed 
only marginally from the standard case (i.e. the case in which the moments are neglected). Thus we found that for 
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FIG. 6: The energy density of the bounce, given by Eq. H57|) . for the W = and W — —0.05 cases as a function of AG'^'^ , with 
the remaining initial conditions being given by Eq. ()77p. The W — case (lower line) is approximately independent of AG'''' 
because at the bounce AG'''' ~ AG^^ , for the system under consideration, however for W (upper line) there is a (small) 
dependence on G'''' . Note that both plots correspond to the same initial range of G'''' = 100 140. The non-vanishing slope 
shows that effects on the bounce density are not just due to the non- vanishing potential. 




FIG. 7: The left hand plot gives the values of the uncertainty conditions, during the evolution, the center one plots the 
corresponding values of AG'''' and Im AG'" . The right hand plot shows Ri, the order of the terms neglected in Eqs. (|2ip - (|22p 
compared to the fractional moments (AG'''')^''^ (the larger values) and (AG^^)^^^ (which happens to be much smaller in this 
evolution) . For this evolution all the remaining moments and correlations remain small (less than « 2%) and the bounce occurs 
at <?ibounco ~ 1.82. The increase of Ri toward the end of the evolution indicates that neglecting the third order moments is no 
longer a good approximation in the region in which the uncertainty conditions (in this case U2) are breached. 



the (large number of) cases investigated, the bounce occurs at approximately the same density, a feature not at all 
unexpected due to the setup of our approximations. The fact that there were small deviations from the standard 
bound is however significant, especially since we show that it is possible for the bounce density to be sometimes 
larger than the standard case. We restricted ourselves to considering only states that remained highly semi-classical 
before, after and during the bounce and still found that the moments noticeably back-react on the trajectories of 
the expectation values. It may be that there are states which are semi-classical at large scales, but which become 
dominated by the evolution of their moments at small scales (i.e. become highly non-semi-classical) and for such states 
the presence of a bounce is not guaranteed by current results in the literature. 
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Appendix A: Evolution equations 

The evolution of the variables V, J and their moments with respect to the scalar field (f> are [40} : 

^ - -:^ + -^±Z_t/3/(i-.)^(^) + 3_^±^^3/(i-.)(^^^^ (Ai) 

6 J + J ^^(2+.)/(i-x)(gyj _ G^J)H/(<^) + 3 2+^ J + J ^^^^ 



l-x{J-J)3 ' ' 2[\-xY{J-Jf 

91/3/(1-2;) o y{2+x)l(\-x) 

d/ / 3 T/(2+2;)/(l-a:) 7 2^3/(1-2:) fT/ i /O^ 



3 v(2+^)/(i-^)j 6F3/(i-^)(^_^;j/2) 



\-X (J - J)3 (J _ J)4 



X\-xf [J-jy \~x {J~Jf 

3 /(2+a;)(l + 2a;)y3^/(i-^)j 2(2 + x) + ft/2) 



G^^W{(j)) 



2 V (l-x)3 J- J ' [l-xf {J- J) 



-(G^'^ -G'^'^ -2G^^)Ty ((/)) + 3 ," = C'Wicj)) 



l-x [J-Jf " ' {\-xf J- J 

4]/3/(l-a:) 



-(G^'^ + G^"'^) + =r^(G^'^ + G^^)W^(0) 

(J - J)^ 



G^^W^(0) - — ' W{(b) 



1-x {J-Jf 2{l-x) (J -J) 

^ ^-G'"^W-(0) - ^/ J^^-^W{cl>) (A3) 



1-X {J-Jf 2(1 -x) (J- J)2 

O T/(2 + 2;)/(l-x) _ _ _ _ 

-((J + J)G^^ - (2y + fi)(G^-' + G^-^) - JG-^^ - JG'''^ 



1-x {J-Jf 
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2 J-J \l-x 3 (J -J) 

+ - ^°")»'(« - 2(1-.) (.-.r ' -^w (*^' 

j^yj 1 _ T/-3/(l-a;) _ o T;-(2+a;)/(l-a:) 

Cicp A ((7 — (7) ± X u — t/ 

+ ^ ""^y ) (_2( J + J)G-^J - 2JG^"^" + 2(V + h/2)(4G^J - 2G^^) + 2JG^^ 
[J-Jf 

A. X i'J ) 



+3(13^ jzrj - -,h J) 



+ 



y3/(l-x) 

-2G^''+4— =-^^'' 

(J - J)2 



6 y(2+3:^)/(l-^) 



G^'^W^(<^) + G'^'^W^(,/>) 

1 — a; J — J 



(A6) 



(J- J)3 
3 i[/(2+x)/(l-x)^(^-) 



(8(y + ?l/2)(2G^^ + G-^-^ - G''-^) - 8JG^-^ - 8JG^-^ + 4?i2y + 2%^) 

6(2 + a;) y(i+2x)/(i-x) 



1 - a; {J-Jf 



{2JG'' - 8{V + h/2)G''-' + 2JG'' + h^J) + 



{l-xf 



J-J 



GJJ . 



Equations of motion for real and iinagiuary parts of J, G^^ and G'^'^ are: 



^(2+.)/(l-.) _ ) 



(ImJ)2 



+ 



3 ]/(4-a;)/(l-a;) 

'i (ImJ)4 
3 



(1-x) 



(ReG-^-^ - G^^")VP^(0) 

y3/(l-a;)^ 

(ImJ)3 



IniG^''l^(^!.) 



3 /(2+^)(l + 2x) 3,/(i_,) 
2^ 2(l-a;)3 



2(im J)2 



y3/(l-a;) 



-1 

^ 4 (Im jf 



4-x y(2+^)/(i-^) 
2(l-a;)2 (ImJ)2 

3 _ 



G*'''H^((/.) 



Re G'' 



+ — 4— -Im G^^y-i - l^^G^^y-^ 
(1 - a;)Im J 4(1 - a;)^ 



(A7) 



(A8) 



d^''"^ = -2(1^ 



y(2+x)/(i-:r)p^gj 3 (2 + x)(l + 2a;) y3a;/(i-a;)ReJ^^^ 



ImJ 



-W^(<^)--- 



(1 - a;)3 



ImJ 



G^'^'Wict)) 



(A9) 
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'2 (Im Jf 
12Im jy-i 



2G 



G 



J J 



Re G 



1 -a: 



(^2Ini G^'''^ 



-Im J 



(A16) 



Appendix B: Fiducial volume 



For the mathematical formulation of cosmological models with non-compact spatial slices, a fiducial volume Vq is 
required in order for the symplectic structure of the form J d^xp^cj) to be defined (i.e. for spatial integrals to be finite). 
Physical quantities must be independent of the value of Vq and in particular remain well defined as Vo — >■ cxd. This 
has, in the past, been the cause of some confusion, so for clarity in this section we describe in detail the Vq dependence 
of the all the parameters in the theory. 

The canonical pair we are using is. 



Uip) c,V} 



3 



8^G7(l-x)/o, 
3 



87r7G ' 8TT-fG 

where the over-tilde implies that the quantity does not depend on Vq. The combination fap^c thus behaves as 



(Bl) 



(B2) 



where /o = foV^. The freedom we have then is to fix the behavior of /o with respect to Vq such that a well-defined 
formulation results. We could, for example, have chosen to set ?/ = 0, as implicitly occurred in some papers on the 
subject. Had we done so, we would have found: 



V = 



3pi 



3pi 



^V 



However, then we have 



8^G7(l-2;)/o 8^G7(l-x)/o ° 



2(l-x)/3 



(B3) 



(B4) 



which only for x = —1/2 is consistent with the primary condition on the scaling behavior of /(p), ensuring that 
holonomies are independent of Vq. In particular, if we were to attempt a limit of large Vq, holonomies or even 
observables such as the density at the bounce would not remain well-defined. 
Instead, we must choose to have /o depend on Vb as. 



(B5) 



so that foP^c = foP^c is independent of the regulator Vq. The parameter fo is Vb-independent only for a; = —1/2, in 
which case fop'^c = foja/a appears in holonomies, as suggested as a special ad-hoc choice by Then we have that 



V = 



3pi 



^Vr 



SttGj (1 - x) fo 



• 



(B6) 



The factor Vb provides the correct dependence on Vb in the Poisson bracket, Eq. (jBl|) . In this case, taking the large- V 
limit corresponds either to taking the large- Vb limit or the small-discreteness limit {fo 0). 

We are interested in physical quantities, such as the kinetic energy density at the bounce (here we focus on the 
W = case, the W case follows similarly), the square-root of which is given by Eq. (15^ . 
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vol 
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Bounce, largeV 
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which is independent of the fiducial vohime for all choices of the lattice refinement parameters x and /q. 
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